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A matched filter maximizes the signal-to- noise ratio of a signal. In the recent work of Corron et al. [Chaos 20, 
023123 (2010)], a matched filter is derived for the chaotic waveforms produced by a piecewise-linear system. 
Motivated by these results, we describe a pseudo- matched filter, which removes noise from the same chaotic 
signal. It consists of a notch filter followed by a first-order, low-pass filter. We compare quantitatively the 
matched filter's performance to that of our pseudo-matched filter using correlation functions in a simulated 
radar appHcation. On average, the pseudo-matched filter performs with a correlation signal-to-noise ratio 
that is 2.0 dB below that of the matched filter. Our pseudo- matched filter, though somewhat inferior in 
comparison to the matched filter, is easily realizable at high speed (> 1 GHz) for potential radar applications. 
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Radar measures the propagation times of trans- 
mitted electromagnetic waves between a target 
and a receiver. A stored copy of the transmit- 
ted waveform is compared to the received sig- 
nal to compute relative distances and locate ob- 
jects. The resolution of a radar is proportional to 
the bandwidth of the transmitted signal and, in 
the case of extremely large bandwidth, the radar 
must be capable of storing and processing large 
amounts of data. One way to maintain high band- 
width and circumvent this sampling limitation is 
to use chaos. Chaotic systems are deterministic, 
non-repeating, and can be engineered to have an 
underlying symbolic dynamics. Rather than sam- 
pling and storing the entire chaotic waveform in a 
radar, one can track the symbols associated with 
the transmitted signal. In this technique, a re- 
ceiver must recover the symbolic dynamics from 
the received signal, which has picked up noise 
from the scattering environment. Recently, Cor- 
ron et al. developed a piecewise-linear chaotic 
system that uses a matched filter to recover the 
system's symbolic dynamics. Using these results, 
we empirically derive a simpler pseudo -matched 
filter that uses high-speed commercially available 
filters to recover the system's symbolic dynamics. 
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I. INTRODUCTION 

A conventional radar system measures the distances 
of targets in the field of view using a signal source, a 
transmitter, and a receiver. In Fig. (TJ a radar trans- 
mitter broadcasts a signal u{t) from the source toward 
an intended target, and the receiver detects a version of 
the transmitted signal that is reflected off of the target. 
Prior to transmission, a copy of the radar's signal is digi- 
tally sampled and stored as Sn- The received signal v{t)^ 
which picks up environmental noise, is filtered and corre- 
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FIG. 1. Example chaos radar. A signal u{t) is transmitted 
to measure target distances and stored digitally as Sn • The 
received signal v{t) is filtered and correlated with Sn- The 
correlation output x(t) peaks at time t = 2rpd + rf, where Tpd 
is the propagation delay from the radar to the target and Tf 
is the delay through the filter. 



lated with Sn- Typical radar signals are non-repeating in 
order to avoid multiple points of correlation and, there- 
fore, the correlation will peak only when the transmitted 
and received signals are aligned. Using the time of the 
output peak in the correlation, the measured range from 
the transmitter to the target is determined. 

The performance of a radar is determined by its abil- 
ity to identify the correlation time between the transmit- 
ted and received signal. In the correlation function, the 
width of the peak scales inversely with the bandwidth of 
the transmitted signal and sets the spatial resolution of 
the radar. In addition, the height of the correlation peak 
above the noise floor, also known as signal-to- noise ratio 
(SNR)^ is determined by the length of the transmitted 
and sampled waveforms as well as the noise from the envi- 
ronment. Thus, the digital storage capacity of the radar 
sets the maximum SNR in the correlation measurement. 
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State-of-the-art radar systems use high-frequency, broad- 
band signals, where the digital sampling and storage of 
the signals can be costly. These radars must balance the 
bandwidth and cost of the systems design while main- 
taining its performance. 

A simple example of an inexpensive, broadband sig- 
nal source is amplified electrical noise. In the past, elec- 
trical noise has been used by radar systems to perform 
ranging measurements ^^i^ The high bandwidth of these 
noise generators yields high-resolution ranging informa- 
tion, but requires fast sampling and large data storage ca- 
pacities. Some recent techniques have been proposed to 
minimize the necessary sampling capacity of noise radars 
using analog delay lines for signal storage But, these 
methods limit the ranging capabilities of the radar. Fur- 
ther research has shown that more efficient radars benefit 
from a structured, rather than a stochastic, signal. 

For instance, various deterministic signal sources have 
been studied in efforts to minimize the necessary data 
storage rate and capacity of a radar. As one exam- 
ple, pseudo-random binary sequences (PRBS's) are often 
used as a radar signal sources. To be implemented as a 
radar signal source, a PRBS is up-converted to a suit- 
able frequency band before transmitting and then down- 
converted at the receiver before correlation.^'^ The main 
advantage of a PRBS is the ability to use one-bit digital 
samplings of the binary sequences, thereby requiring low 
amounts of data-storage capacity. This allows for longer 
sequences of the transmitted waveform to be stored, thus 
enhancing the radar's SNR without increasing the cost 
of the system. The main disadvantage of a PRBS is 
that it requires computational power to generate and its 
sequence eventually repeats, which ultimately limits its 
performance. Many other radar concepts like this one 
exist, each with advantages and disadvantages, and to- 
day the radar community continues to develop broadband 
signal sources.^'- 

One novel approach to a radar is to use chaotic wave- 
form generators, which are believed to have several prop- 
erties that make them ideally suited as signal sources 
for radar applications One defining feature of a chaotic 
system is that it can generate a signal that does not 
repeat in time. Chaotic signals are also often inher- 
ently broadband. High-speed chaos has been observed 
in optical and electronic systems with frequency band- 
widths that stretch across several gigahertz. Such 
broadband chaos has been studied for its applications in 
high-resolution ranging and in imaging — These appli- 
cations use the non-repeating aspects of the high-speed 
chaos. 

However, conventional chaos radars do not take full 
advantage of the chaotic signal source. In addition to 
producing broadband, non-repeating signals, chaotic sys- 
tems are deterministic and extremely sensitive to small 
perturbations. By not using these properties, chaos 
radars add no benefit over noise radars, requiring high- 
sampling to perform correlations. Thus, to take full ad- 
vantage of chaotic systems in a radar application, a chaos 



radar needs to benefit from the determinism or sensitivity 
of the chaos in addition to its noise-like properties. 

Recently, Corron et al proposed a novel chaos radar 
concept that uses the dynamics produced by a piecewise- 
linear harmonic oscillator J^— It produces simultane- 
ously a chaotic signal and a binary switching state that 
completely characterizes its chaotic dynamics. In the 
proposed radar, the chaotic signal serves as the signal 
source (see Fig. [T]) and is transmitted, while a copy of a 
switching state is stored using a one-bit digital sampling. 
For a radar receiver, Corron et al. derived the analytical 
form of a filter that is matched to the chaos produced 
by the system. A matched filter is a linear operation 
that optimizes the SNR of a signal in the presence of 
additive white Gaussian noise (AWGN). Their matched 
filter recovers the switching states from the received sig- 
nal, which are then passed to the correlation operation 
of Fig. [H This technique uses deterministic aspects of 
the chaos, making it an improvement over a noise radar. 
With reduced data storage and an optimal SNR^ their 
architecture could effectively reduce the cost of a radar 
system. 

Corron et al. implement their piecewise-linear design 
using an LRC (inductance-resistance-capacitance) oscil- 
lator that operates in the kHz frequency range. It is 
difficult to a realize high-frequency version (> 1 GHz) 
of this system because of parasitic capacitances, induc- 
tances, and the inherent time delays in the propagation 
of signals in LRC circuits. As it stands, there is no high- 
frequencv realization of the piecewise-linear system from 
Ref. [13 or the associated matched filter. In order to 
have resolutions that are comparable to state-of-the art 
radar systems, the waveforms and switching states pro- 
duced by this chaotic system must be scaled to higher- 
frequencies and to broader bandwidths. Thus, techniques 
for increasing the speed of this system are of interest. 

As a first step toward a high-frequency implementation 
of Corron et al. ^s system, we compare the performance of 
the matched filter presented in Ref. [17] to a set of stan- 
dard filters (first order low-pass filter and notch filter). 
These standard filters operate at high-frequencies, are in- 
expensive, well characterized, and readily available. Cas- 
cading these standard filters allows us to realize a pseudo- 
matched filter for the chaos produced by the piecewise- 
linear harmonic oscillator. We define a pseudo-matched 
filter as a sub-optimal linear operation (when compared 
to the matched filter) that removes AWGN from the 
corresponding (matched) waveform and performs com- 
parably to the matched filter for the system. We are 
interested in a chaos-based radar system such as that de- 
scribed by Ref. [HI that could benefit from integrating 
readily available filters for high-speed applications. 



II. MATCHED FILTER REVIEW 

To motivate our analysis, we briefly review the char- 
acteristics of the chaotic system and matched filter pre- 
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sented in Ref. [17| within the context of a radar apphca- 
tion. Consider a harmonic oscihator with negative damp- 
ing — /3 and with a piece- wise constant driving term s{t) 
whose behavior is governed by the differential equation 

u{t) - 2pu{t) + {ul + p^){u{t) - s{t)) = 0, (1) 

together with a guard condition on the output variable 
u{t) that switches the sign of s{t) 
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)>o 
)<o 
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Figure [2] shows a time series of the variables u{t) and 
s{t) with the parameters cJq = 27r and /3 = ln{2). We in- 
tegrate Eqs. (P)-© using MATLAB's ODE45, where the 
switching condition is monitored using the integrator's 
event-detection algorithm. The attractor for this system 
is low dimensional and is plotted in Fig.[2]3. The dynam- 
ics can also be viewed as a chaotic shift map^^-^ as seen 
in Fig.[2t. For this system, u{t) oscillates with a growing 
amplitude and fixed oscillation frequency /o = CjOq /27r 
about a piece- wise constant line s{t). The switching 
times of s{t) depend on the local maxima and minima 
of u{t), and the times between the local maxima and 
minima of u{t) are fixed by the fundamental frequency 
/o- Thus, the maximum rate of switching in s{t) is lim- 
ited to l//o. Using a one-bit digital sampling of s{t) at 
a sampling frequency that is greater than or equal to /o, 
we are able to store a record of all switching state values 
Sn- Similar to the case for a PSRB, information about 
the transmitted waveform can be stored with minimal 
sampling and memory, enhancing the potential SNR of 
a radar correlation measurement. 

In addition to its data storage capabilities, chaos from 
this system can be further exploited using a matched 
filter. Corron et at. demonstrate that the switching in- 
formation s{t) can be recovered using a matched filter at 
the radar receiver. The match filter is described by 



y{t)=v{t-2^/ujo)-v{t), 



(3) 



(mit) + 2pimit) + {ujI + P^)Uit) = {ujI + P^)y{t), (4) 

where v{t) is the input signal and ^m(^) is the analog out- 
put of the matched filter. In Fig. [3l we examine the out- 
put of the matched filter when it is driven by v{t) = u{t) 
+ AWGN. The original signals u{t) and s{t) are plotted 
in Figs. [3K-b. The switching state s{t) is plotted with a 
one-bit digital sampling Sn at a sampling frequency /q. In 
Fig. [St, the time series v{t) has a SNR of —5.9 dB, where 

SNR = 101ogio(^A/'i^input), ^A^i^input = ^l/^AWGN^ 

and where and ct^^^qj^ are the input signal u{t) and 
additive noise powers, respectively (see the Appendix for 
details on the additive noise). The matched filter is op- 
timized to remove the noise from v{t) and recover the 
binary switching information Sn for the correlation oper- 
ation. 
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FIG. 2. Chaos from a piecewise- linear harmonic oscillator 
with negative damping, (a) Time series of the analog vari- 
able u{t) (green) and the nonlinear switching state s{t) (blue 
dashed line), (b) Chaotic attractor in phase space, (c) 
Chaotic shift map created by sampling u{t) using the times 
t* from Eq. (2), where Vn = ^^(t*) if |ix(t*) - 1| < 0. 



We compare Sn to the solution for Eqs. (|3]) - dH, ^m(^), 
shown in Fig. [3]i. The solution ^m(^) is the output of 
the matched filter and is a nearly noise-free waveform 
that transitions approximately between two states, one 
positive and negative (defined by a dotted black fine at 
^m(^) =0)- In between these transitions, ^m(^) oscil- 
lates with a small, irregular amplitude at a frequency 
approximately equal to /q. Upon further examination, 
the matched filter output signal £,m{t) follows the same 
transitions from the original switching state s{t). Us- 
ing a correlation between s{t) and ^m(^), we determine 
the time delay through the matched filter to be ap- 
proximately l/(2/o). We compensate for the delay, and 
sample ^m(^) at /o, assigning binary values using the re- 
lation: — 1 if ^m(^n) ^ and +1 if ^m(^n) > 0, where tn 
is the n*^ sampling time. These binary values, shown in 
Fig.[3]i, match s^, demonstrating that the matched filter 
has recovered the switching information from v{t). 

In a simulated radar application, Corron et ai use a 
tapped delay line to perform a time-domain correlation 
operation between the digitally stored Sn and the analog 
output of the matched filter ^m(^)- The tapped delay line 
is described by 



N 



(5) 



n=l 



where A^ is the length of the stored sequence and Xm(^) 
is the output of the correlator (we have used the linear- 
ity of Eqs. (d]) - (jlj) to rearrange the operations in Ref. 
[l8l|). The tapped delay line is a convenient technique for 
calculating in real time the correlation between a binary 
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FIG. 4. Schematic representation of the tapped delay hne in 
Eq. dSJ. The switching state s{t) is sampled and Sn is stored 
for n = 1 to n = A^, where, N = 100 in this case. The output 
of the matched filter ^m(t) drives the tapped delay line and 
the output sum Xm(t) peaks at time t — 2Tpd + Tm- 



receiver. The matched filter for chaos, combined with 
a tapped delay line, provides an architecture for rela- 
tively quick and inexpensive correlations between a bi- 
nary switching signal and a recovered analog signal, both 
generated from a single chaotic system. These benefits 
present a platform for a high-performance radar with im- 
proved SNR over conventional chaos radars. 



FIG. 3. Temporal evolution of (a) n(t), (b) s(t), (c) v{t) with 
SNR = -5.9 dB, and (d) ^m(t). The horizontal axes in (b) 
and (c) are shifted by Tm and (2Tpd + Tm), respectively, where 
Tpd is the propagation distance to an intended target and 
Tm is the time delay through the matched filter. Above the 
signals ^m(t) and s(t) (blue), a single-bit discrete sampling of 
the waveforms (red dots) is shown. 



and an analog signal. 

To better understand the operations performed by a 
tapped delay line, we use a pictorial representation of Eq. 
([5]). In Fig.m the output of the matched filter ^m(^) splits 
into N copies, each of which are successively delayed by 
Tn/o = ^, where n is an integer. The resulting copies are 
multiplied by the corresponding stored Sn and summed 
continuously in time. The output Xm(^) peaks 2rpd + Tm, 
where Tpd and Tm are the propagation delays of the signal 
to the target and through the matched filter, respectively. 
In the context of Fig. (TJ the correlation operation of this 
chaos radar can be performed between a transmitted and 
received signal using the tapped delay line in Eq. (j5j). 
Thus, using the chaotic waveform generator from Ref. 
f\j\ , combined with the matched filter and tapped delay 
line, we arrive at a chaos radar. 

This particular chaos radar benefits in two ways from 
the deterministic characteristics of the system's chaos. 
The first benefit is the link between the non-repeating 
waveform u{t) and the switching state s{t): A binary 
sampling of s{t) can completely characterize the dynam- 
ics in u{t). A second benefit is the ability to derive a 
matched filter that optimizes the output SNR of the 



III. MATCHED FILTER ANALYSIS 

With this background and motivation, we now examine 
Corron et a/.'s matched filter for chaos in the frequency 
domain. Using the transfer functions of Eqs. ([I]), (j3j), and 
(jH), we examine the spectral properties of the matched 
filter. For the purposes of our analysis, we focus on the 
magnitudes of transfer functions. Because the phase of 
the matched filter is approximately linear with frequency 
/ when / < /o, it preserves the transition information 
in Sn- We will derive empirically a pseudo-matched filter 
using a combination of standard filters and that also pre- 
serves the transition information in Sn- Lastly, we com- 
pare the pseudo-matched filter performance to the true 
matched filter for chaos in a simulated radar application. 

First, we analyze the spectral properties of the chaotic 
dynamics from u{t) and s{t) as well as the driving sig- 
nal v{t). In Fig. [51 we plot the spectral amplitudes 
for u{t), s{t), and v{t), where v{t) = u{t) + AWGN 
{SNR = -5.9 dB). One should take note that, in Fig. [5^, 
there is no maxima in the frequency spectrum at the 
fundamental frequency of oscillation, because the phase 
of u{t) switches by tt each time s{t) switches states. This 
demonstrates that, if the spectrum of u{t) is scaled- up 
in high-frequency system (> 1 GHz), the bandwidth in 
u{t) would stretch over several gigahertz. In addition to 
its radar properties, a high-frequency broadband spec- 
trum in u{t) would provide a useful carrier signal for low- 
profile or ultra- wideband technologies In Fig. [St, the 
frequency spectrum of the AWGN in v{t) covers infor- 
mation about u{t). The matched filter is engineered to 
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FIG. 5. Frequency spectra for the chaotic time series of (a) 
u{t), (b) s(t), and (c) v{t). The frequency spectra are cal- 
culated using s, and v, the discrete Fourier transforms of 
u(t), s(t), and v{t), respectively. In these plots, the frequency 
axes have a spectral resolution of 10~^, and the spectral am- 
plitudes have been averaged over a window of 10~^. 



process the spectrum of v{t) to recover Sn- 

Next, we examine the transfer function of the matched 
filter. We take the Fourier transform of Eqs. (|3|) - 
and obtain the transfer functions i^input and Hq. We 
combine these transfer functions to obtain the transfer 
function for the matched filter H^n- They read 



^inputi^j — ;r = - 



1 



27Tiiy 



47r2 + /32 



-^input-^07 



(6) 



(7) 



(8) 



where v = j j jo and ^m, ^, and v are the Fourier trans- 
forms of ^m(^), v{^\ and v(t\ respectively. We plot the 
magnitudes of i^input, and as a function of fre- 
quency V in Fig. [6l In Fig. [6t, we also plot the phase 
of i^m, where the phase is approximately linear with fre- 
quency for < 1 and thus preserves timing information 
in Sn. 

We analyze i^input and individually to better un- 
derstand the matched filter's transfer function. We 
factorize i^input into two linear operations (i^^input = 
^notch ^integrator), a uotch filter and an integrator 



i^notch(z^) = - = 



(9) 



OQ 

■o 



OQ 

■o 



OQ 

E 
X 



■o 




FIG. 6. Magnitudes the transfer functions (a) i^input, (b) ifo, 
and (c) ifm. In (d), the phase of ifm is given as a function 
frequency. 



-^integrator v^j — ~ — ^ . i 



(10) 



where q is the Fourier transform of an intermediate input- 
output variable and i^notch and i^integrator are the trans- 
fer functions of a notch filter and integrator, respec- 
tively. We plot the magnitudes of i^notch and i^integrator 
in Fig. [3 In Fig. [7^, the magnitude of the notch fil- 
ter's transfer function goes to minus infinity at integer 
multiples of the fundamental frequency /q. Similar types 
of filters have been used previously in chaotic systems 
to stabilize periodic orbits using continuous-time control 
methods*21i2^ The transfer function in Eq. (p!Q|) is a stan- 
dard operation for integration. In Fig. [T^, the function 
I ^integrator I divcrgcs to infinity dit u = and falls off at 
a rate of l/v with increasing v. When cascaded, i^notch 
and i^integrator Complement one another to form a filter 
that preserves low frequencies, eliminates integer multi- 
ples of /o, and cuts out high-frequencies (see Fig. [6ti). 

Upon inspection of Rq in Eq. ([7j), we see that it takes 
on a functional form that is equivalent to the Fourier 
transform of the dynamical system in Eq. ([T]), given by 
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FIG. 7. Magnitudes of the transfer functions (a) i/notch and 



(b) H, 



integrator • 



where u and s are the Fourier transforms of u{t) and 
respectively. Thus, Eq. ([7]) contains specific spec- 
tral information about the system's dynamics. When the 
transfer function is applied with i^input? it reshapes 
the spectrum of the transfer function for the matched fil- 
ter i^m near = 1, as seen best in Fig. [6t. Thus, the 
operations of the matched filter for chaos can be sepa- 
rated into four criteria: (i) a notch filter that eliminates 
the fundamental oscillation frequency /o, (ii) an integra- 
tor that preserves the low frequencies and (iii) cuts off 
high frequencies falling off as 1/z^, and (iv) a dynamical 
filter that reshapes the transfer function near the funda- 
mental oscillation frequency /q. These four criteria are 
the foundation for our derivation of the pseudo-matched 
filter for chaos. 



IV. PSEUDO-MATCHED FILTER 



^shifted-notch(^) = = "(1 + e""*^), (12) 

where Vin and Vout are the Fourier transforms of the in- 
put signal Vin and output signal Vout^ respectively. We 
plot the magnitude of i^shifted-notch in Fig. ^ (compare 
to i^notch from Fig. [TH). In both plots, the fundamental 
frequency /o is blocked. However, in Fig. [8^, the lower 
frequencies {v < 0.5) are not cut. Thus, the shifted- 
notch filter performs two of the four operations from the 
matched filter; (i) it eliminates the fundamental oscilla- 
tion frequency fo and (ii) preserves low frequencies. 




Our strategy for designing a pseudo-matched filter for 
chaos is to simplify the four criteria of the matched filter 
using transfer functions from components that are readily 
available at high-speed. In the upcoming section, we sat- 
isfy criteria (i) and (ii) using a single transfer function. 
We also show that criterion (iii) is more easily accom- 
plished without an integrator. Lastly, we demonstrate 
that criterion (iv) is not necessary for our appHcations. 

To begin constructing our pseudo-matched filter, we 
select a different notch filter that is shifted in frequency 
but still blocks the fundamental frequency /q. Most 
notch filters block integer multiples (n = 0, 1, 2, 3, ...) of a 
single frequency. Instead, we choose a notch filter that is 
shifted to block odd integer multiples (2n + l = 1, 3, 5, ...) 
of a single frequency. Since the matched filter attenuates 
frequencies above v = 1, we conjecture that the only im- 
portant spectral notch is at /o, and all high-order even 
notches are not included in our pseudo-matched filter. 
The transfer function of our shifted-notch filter is 



FIG. 8. Magnitudes of the transfer functions (a) i^shifted-notch , 
(b) i/iow-pass, (c) Hp. In (d), the phase of Hp is given as a 
function frequency. 

In the time domain, the shifted-notch filter of Eq. ([12]) 
is expressed by 

^out(^) = ^(^in(^) + ^in(^ - Tt/cJo)), (13) 

We compare Eq. ([13]) to Eq. ([3]) and note that the out- 
put is no longer related to the input through a derivative. 
Also, the time-shift on the input signal is halved (tt/cJo 
instead of 27r/cJo) and the shifted input Vin{t — tt/uJo) is 
summed with the present state Vin{t). In an experimental 
setting using high-speed electronics, where and Vo^t 
are voltages, this shifted-notch filter can be realized using 
a voltage divider, time-delays (realized, for example, by 
coaxial cables), and an isolating hybrid junction, as illus- 
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trated in Fig.[9l The lengths of the two cables used in this 
realization of the filter are chosen such that the difference 
in propagation times for electromagnetic waves to prop- 
agate through them is: tb — ta = tt/cJo- The isolating 
hybrid junction sums the outputs: v{t — ta) -\- v{t — tb)- 
We shift time t ^ t -\- ta to arrive at the output signal 
Vout in Eq. (14). This realization of the shifted- notch 
filter can scale to high-speed voltages (> 1 GHz). 
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(15) 



where Vin and Vout are the Fourier transforms of the in- 
put and output signals of the filter, respectively. We plot 
the magnitude and phase of in Fig. [Hb. In the fig- 
ure, the phase of is approximately linear and thus 
preserves timing information from v(t). For comparison 
to the matched filter, see Fig. [6b. Qualitatively, the two 
filters follow similar trends in both magnitude and phase. 
The phase in the pseudo-matched filter has a lower slope 
in its frequency dependence; a lower slope just consti- 
tutes a shorter time delay through the filter. However, it 
is clear by comparison of the magnitudes that the pseudo- 
matched filter is not performing the same operations as 
the matched filter. 



FIG. 9. Pictorial realization for a high-speed (> 1 GHz) 
shifted-notch filter for voltae: es V'lYi and ^;out. An example of a 
broadband, high-frequency power splitter is the Mini-Circuits 
ZFRSC-42-S, and an example of a broadband, high-frequency 
hybrid-junction is the M/A-COM H-9. 



Continuing the construction of the pseudo-matched fil- 
ter, we use a first-order low-pass filter to attenuate high 
frequencies, rather than an integrator. We avoid the need 
for an integrator because the shifted-notch does not cut 
off low frequencies. The transfer function of the low-pass 
filter is 



ow-pass 



^out 



1 + 2mv ' 



(14) 



where xin and Xout are the Fourier transforms of the in- 
put and output signals, respectively, and the low-pass 
cutoff frequency is = /l//o- We plot the magnitude 
of i^iow-pass in Fig. [Hb. In the figure, i^iow-pass leaves the 
spectral amplitude of frequencies below vl unchanged, 
while suppressing frequencies above z^l- Beyond z^ = z^l, 



the rate of the spectral roll-off of |i^io 



is not V , 



but rather ^ (1 -\- . A first-order low-pass filter is 
a standard electronic component for filtering an input 
voltage to obtain an output voltage Xout and satis- 
fies approximately the third component of the matched 
filter criteria (iii). We note that higher-order low-pass 
filters (Butterworth, Chebyshev, etc.) are also available 
at high-speed. 

When constructing our pseudo-matched filter, we ne- 
glect the dynamical filter that reshapes the spectrum (iv). 
We show that using just the shifted-notch and low-pass 
filters allows us to achieve comparable performance to 
the true matched filter in a simulated radar application. 
Thus, we cascade the shifted-notch and low-pass filters 
(i^shifted-notch ^low-pass) to arrive at the transfer function 



We now apply the pseudo-matched filter to the chaotic 
waveform generated by Eqs. ([I])-© and examine its out- 
put in the time-domain. We drive the pseudo-matched 
filter with v{t) = u{t) + AWGN, where v{t) has a SNR 
of —5.9 dB (see Fig. [3b). In Fig.fTOk. we plot the pseudo- 
matched filter's output ^p(t). From the figure, we see that 
pseudo-matched filter has effectively removed the main 
oscillation frequency /o and what remains is a digital-like 
signal ^p(t), where the noise has also been reduced. We 
note that a considerable amount of noise is still present 
in comparison to Fig. [Sji. Using a correlation between 
the original s{t) and ^p(t), we determine the time delay 
through the pseudo-matched filter Tp to be approximately 
O.I4//0. We compensate for the delay, and sample ^p(t) 
at /o, assigning binary values using the relation: —1 if 
^p(^n) ^ and +1 if ^p(tn) > 0, where in is the n*^ sam- 
pling time. In Fig. [TOk . we see that, with this particular 
SNR, the discrete samphng of ^p(t) is equivalent to Sn 
from Fig. [8)3. 



Next, we parallel the construction of the time delay 
tap from Eq. ([5|) for the output of the pseudo-matched 
filter. In this case, the time delay tap is 



N 



(16) 



n=l 



where Xp(^) is the output of the time-domain correla- 
tion. We plot an example of Xp(^) in Fig. [TOb using the 
the same chaos and Sn that were used to calculate Xm(^) 
in Fig. [4)3. By visually comparing Xp(^) to Xm(^), we see 
that the output correlation peaks are qualitatively simi- 
lar, but Xp(^) has more noise. In the remaining section, 
we establish criteria for quantitatively comparing these 
correlation waveforms and use these criteria to weigh 
each filter's performance. 
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FIG. 10. Output of pseudo-matched filter, (a) Time series of 
the output of the matched filter ^p(t) (blue) while driven by 
v{t) = u{t) + AWGN {SNR = -5.9 dB). The signal ^p(t) is 
sampled with uniform spacing (red dots) at a clock frequency 
/o = cJo/Stt. Above the waveform, a single-bit discrete sam- 
pling of the waveform is shown. From the figure, we see that 
all of the relevant information from Sn is encoded in ^p(t). 
(b) The switching state s(t) is sampled and Sn is stored for 
N = 100. The output of the pseudo-matched filter ^p(t) drives 
Eq. fT6|) and the output Xp(^) peaks at time t = 2Tpd + Tp. 



Next, we measure the output SNR^s of the matched 
and pseudo-matched filter correlations using the correla- 
tion peak heights am,p and the surrounding correlation 
noise floors. The output SNR in Xm,p(^) is 



SNRr, 



m,p 



(17) 



N|m,p 



where am,p is the peak height of ^m,p(^) from a Gaus- 
sian fit (see Appendix) and cr^|^ p is the output variance 
of the correlation noise floor (note that the mean of the 
noise floor ~ 0) for the matched and pseudo-matched fil- 
ters, respectively. We present a summary of these quan- 
tities in the block diagram shown in Fig. [TTh . In the di- 
agram, we also review the waveforms and processes used 
for generating ^m,p(^) and Xm,p(^) and highlight the two 
relevant quantities, SNRi^put and SNRm^p- We calcu- 
late SNR 

m.vy as a function of the input SNRinput- The 
results of these calculations are given in Fig. [TTb . In ad- 
dition, we use the distributions from ^m(^) and ^p(t) from 
the two different cases v{t) = AWGN and v{t) = u{t) to 
derive a an analytical prediction for SNR^^p (see Ap- 
pendix for derivations). These theoretical predictions 
are plotted with dotted lines in Fig. [TTb . These plots 
represent the performances of the matched and pseudo- 
matched filters in a simulated radar. 



V. MATCHED VS. PSEUDO-MATCHED 

In radar applications, the ability to correctly identify 
the location of the correlation peak in the correlation 
operation is the useful measure. Therefore, in order to 
compare quantitatively the output from the matched fil- 
ter and pseudo-matched filter, we weigh each filter's per- 
formance based on the peak width and output SNR of 
Xm,p(^)- We also present an approximate analytical form 
for each correlation's output SNR. 

The peak widths of the output-correlation functions 
give the resolutions of each radar system. We measure 
Am and Ap, the full- width at half maximum (FWHM) 
time of the correlation output peaks using the matched 
and pseudo-matched filters, respectively. For the most 
ideal measure of each filter's correlation peak width, we 
measure Am,p in cases where no noise is present in the re- 
ceived waveform v{t) = u{t). We note that these widths 
are independent of TV, the number of stored data points 
in the correlation calculation of Eqs. (|5]) and (fT6|) . Us- 
ing a Gaussian fit to the peak of Xm,p(^), we obtain peak 
widths Am/o = 0.55 and Ap/o = 0.73 (see Appendix). 
Using these values of Ani,p and scaling /o to 1 GHz, we 
calculate the theoretical resolutions of the matched and 
pseudo-matched filters to be 0.17 m and 0.22 m, respec- 
tively. In this example, the ranging resolutions differ by 
5 cm. Thus, this is not a critical difference for radar ap- 
plications that localize targets like planes or cars, and 
the pseudo-matched filter has an acceptable ranging res- 
olution in comparison to the matched filter. 



From Fig. [TTb . it is clear that the matched filter out- 
performs the pseudo-matched filter in the output SNR 
of a radar correlation. Without noise in the system, 
the matched and pseudo-matched filters perform with 
output correlation SNR's of 2.6 + lOlogio(A^) dB and 
1.3+101ogio(iV) dB, respectively. For SNRi^put = 1/100, 
the output SNRs decrease to -2.0 + lOlogio(A^) dB and 
— 4.4 + 101ogxo(^) dB, respectively. In Fig. [TTb. the aver- 
age difference between SNR^ and SNRp is 2.0 dB. We 
note that this difference is independent of N and there- 
fore fully characterizes the filter performances. Thus, 
where small loss is acceptable in the performance of the 
radar, the pseudo-matched filter is a simpler alternative 
to the system's analytically matched filter for chaos. 

As a final example, we use the theoretical SNR^^p 
to predict when the matched and pseudo-matched fil- 
ters will fail in a radar application. Failure occurs when 
SNR^^p falls below a certain threshold. For a radar sys- 
tem that is capable of storing TV = 50 data points and has 
a desired output correlation SNR of 33 dB, the matched 
and pseudo-matched filters will fail at a 1 /S* A/^i^input 

of 

approximately 25 and 70, respectively. If, in this applica- 
tion, the input SNR is such that 1/ SN Ri^itput < 10, then 
a radar with either the matched or pseudo-matched filter 
will be able to range, on average, without failure. The 
choice between the matched and pseudo-matched filter is 
therefore an application-dependent problem, and, as the 
bandwidth of this system scales higher, one must also 
begin to weigh each filter's high-frequency capabilities as 
well as its baseline performance. 
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FIG. 11. (a) Block diagram for testing the matched and 
pseudo-matched filters in a simulated radar application, (b) 
Output-correlation SNR^s of the matched (blue □) and 
pseudo-matched (red Q) filters scaled by A/" on a logarith- 
mic scale as a function of 1/5* A/'i?input • For each value of 
SN Rinput, 100 calculations of SNRm,p were performed using 
sequence of Sn for n = rio to n — Uq -\- where rio is a 
random positive integer and = 50. The mean value of the 
calculated SNRm,p is plotted with the respective standard 
deviations. The bue and red dotted lines give the theoretical 
predictions of the SNRm,p as references for the matched and 
pseudo-matched filters, respectively. Cases for larger N were 
verified to have quantitatively similar results. 



VI. CONCLUSIONS 

In conclusion, for the chaotic system presented in Ref. 
[13 , we derive empirically a sub-optimal filter for remov- 
ing noise from the waveforms generated by this dynam- 
ical system. This sub-optimal filter performs approxi- 
mately three out of four of the linear operations from the 
matched filter: (i) eliminates the fundamental oscillation 
frequency /o, (ii) preserves the low frequencies, and (iii) 
cuts off high frequencies. Our filter, deemed a pseudo- 
matched filter, is composed of a shifted-notch filter and a 
first-order low-pass filter. In the context of a radar con- 
cept that uses a time delay tap as a correlation measure, 
we have shown that, depending on the application, the 
pseudo-matched filter may be an acceptable substitute 
for the matched filter. In addition, we acknowledge that 
the psuedo-matched filter can be further improved using 
higher order low-pass filters and additional shifted-notch 
filters. We present this current version of the pseudo- 
matched filter to illustrate our method for its derivation 



and emphasize its simplicity. 

Lastly, we note that our analysis highlights the fiex- 
ibility of Corron et a/.'s findings. The chaos from the 
dynamical system in Eqs. ([I]) - (|2l) can be processed 
by a linear filter to recover an underlying digital wave- 
form. We capitalize on this system's elegance to create a 
pseudo- matched filter. Although it is less optimal when 
compared to the matched filter, the pseudo-matched fil- 
ter shows that the advantages of this chaotic system can 
be adapted for applied settings that use commercially 
available, high-speed filters. 
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Appendix A: Additive Noise 

Because the output from MATLAB's ODE45 uses a 
variable timestep, we resample u{t) using a linear inter- 
polation with time steps St^ where Stfo = 10~^ . To sim- 
ulate environmental noise, we add noise to the waveform 
u{t) using random numbers spaced by time units St- The 
random numbers are calculated from a Gaussian distri- 
bution with zero mean. For different points along the 
l/SNRinput axis of Fig.fTTb. the variance of the AWGN 
is varied accordingly. 



Appendix B: Gaussian Fits 

We measure the correlation peak width and height us- 
ing a Gaussian fit 

/(t)=a^,pe(*-(2-p<^+^'".p))^/2<p, (Bl) 

where am,p and Cm,p are free parameters that are fit to 
the correlation peak heights and widths. Using f{t) to 
fit Xm,p(^), we obtain a FWHM peak width Am,p/o = 
2^/2ln{2)cm,p and peak height am, p. 



Appendix C: Analytical SNR's 

We derive analytical forms for the output-correlation 
SNR of the matched and pseudo-matched filters. To do 
so, we approximate Eq. (fT7|) as 

N m,p 1 m,P 2 m,p 
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where Am,p is a constant that characterizes the growth 
rate of the correlation peak height with a^^^ p ^ 
constant determined in the noise- free case where v{t) = 
u{t), and cr2|mp ^ function of SNRinput in the case 
where v{t) = AWGN. Recall that the numerators and 
denominators of Eq. (|Cip represent the power of peak 
heights of the correlation and the surrounding noise floor, 
respectively. We derive each of the three terms Am,p, 
^i|m p' ^2|m p following scctious. 

The correlation peak heights for the matched and 
pseudo-matched filters grow at different rates. In the cor- 
relation operations of Eqs. (|5|) and ([16]), the peaks occur 
at times p = 2rpd + rm,p for the matched and pseudo- 
matched filters, respectively. At time t^ p, ^m,p(^) is 
aligned with Sn and the output correlation is 

N 

Xm,p(C,p) ~ E l?m-p(C,p - ^n) I ~ ^m,piV. (C2) 

n=l 

We approximate Am,p from the local maxima of |^m,p(^)|- 
To do so, we examine the noise- free case where v{t) = 

u{t) and collect a subset of points |^m,p(^m,p)| where tm,p 
are the times of local maxima in |^m,p(^)|- We average 

|^m,p(^m,p)| to obtain Am = 0.67 and = 0.51 using a 
time- length tfo ~ 10^. 

To approximate the value of o-^^^^, we also exam- 
ine ^m,p(^) in the noise- free case where v{t) = u(t). 
The deterministic noise floor in a correlation measure- 
ment is also known as its side- lobes; the side- lobes re- 
sult from non-zero contributing terms in the correlation 
Xm,p(^) when t ^ ^m,p- Using the central limit theo- 
rem, we approximate the variance of these nonzero terms 
as (7^1^ pA/", where ^ is the variance of the signal 

^m,p(^)- In this approximation, we find that cr^^ = 0.25 
and cr^ip = 0.20. 

It remains to calculate the contributions to the noise 
floor of the correlation from additive noise. To do so, we 
examine Cm,p(^) in the case where v{t) = AWGN. Simi- 
lar to the case for the side- lobes, we use the central limit 
theorem to approximate the contribution of the AWGN 
to the noise floor of the correlation as cr2\m p^' where 
^2|mp variance of the signal ^m,p(^)- However, 

^2|m p depends on the variance of the input AWGN 

^2|m,p = ^m^pO-^-^GN^ (^3) 

where am,p is the noise attenuation factor of the matched 
and pseudo- matched filters, respectively. We measure the 
values am = 1/289 and = 1/255. Lastly, we use that 



^AWGN = ^l/{SNRinpnt) to rewrite Eq. ([Cl]) as 

SNRm,^ ^ N- (C4) 

^l|m,p + ^m,p5^^^p^^ 

where = 1.34 is the power of the chaotic signal u{t). 
We plot Eq. ()C4p as a function of l/SNRmput for the 
matched and pseudo-matched filters in Fig. [TTb . 
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